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■ Abstract 

It is assumed that during inflation, all energy was contained in a slow- 
rolling inflaton field (p. The particles constituting the Universe are created due 
•/^ . to interactions with the field coherently oscillating after infiation. The lead- 

^ . ing channel of the particle production would be the non-perturbative regime 

^ . of parametric resonance, preheating. Some of the recent developments in the 

00 ' theory of preheating are briefiy reviewed. As a prototype we use the model 

Q> A^4 _|_ 1^02 _|_ expanding universe. In different domains of parame- 
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ters the character of preheating (0 x) (0 ~^ ^'P) is different, ranging from 
regular Lame to stochastic parametric resonances. Fortunately, in many impor- 



D . tant cases simple analytic resonant solutions are found. Eventually, a picture 

which unifies the different regimes emerges. 

^ ■ 1. Creation of Particles from Inflatons 



In the inflationary scenario the Universe initially expands quasi-exponentially 
in a vacuum-like state with vanishing entropy and particle number densities. At 
the stage of inflation, all energy is concentrated in a classical slowly moving in- 
flaton fleld (j). Consider a simple chaotic inflation. The fundamental Lagrangian 
£(0, X, '0, Ai, hik, ...) contains the inflaton part with the potential V{(f)) and other 
flelds (0, X, "0, ^j, ^jfc, •••) which give subdominant contributions to gravity. The 
Friedmann equation for a scalar factor ait): H^ = 3^^|0^ + ^(0)^ (with 

H = a /a) and the Klein-Gordon equation for 0(t): + 3if0 + V^^ = deflne 
the evolution of the background flelds. Soon after the end of inflation, an almost 
homogeneous inflaton fleld 0(t) coherently oscillates with a very large amplitude 
of the order of the Planck mass ~ O.lMp. This scalar fleld can be considered 
as a coherent superposition of 0-quasiparticles with zero momenta, inflatons at 
rest. The amplitude of oscillations gradually decreases not only because of the 
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expansion of the universe, but also because of the energy transfer to particles 
created by the oscillating field. At this stage we shall invoke the rest of the 
fundamental Lagrangian C{(j),x,4',Ai,hik,...) which includes the other bosons 
and fermions and their interactions with the inflaton field, as well as the self- 
interaction of the inflaton field. These interactions lead to the creation of many 
ultra-relativistic particles from infiatons. In this scenario all the matter consti- 
tuting the universe would be created from this process of reheating. A toy model 
describing the interaction between infiatons and other massless Bose particles x 
will be £ = —^g'^(p'^x'^- Gradually, the inflaton field decays and transfers all of 
its energy non-adiabatically to the created particles. If the creation of particles 
is sufficiently slow, the particles would simultaneously interact with each other 
and come to a state of thermal equilibrium at the reheating temperature T^,. This 
gradual reheating can be treated with the perturbative theory of particle creation 
and thermalization. That elementary theory of reheating was developed long 
ago with the first models of infiation. Consider the Heisenberg representation 
of the quantum scalar field x, with the eigenfunctions Xk{t) e~*^^ with comoving 
momenta k. The temporal part of the eigenfunction obeys the equation 



Xk + 5-Xk +[^-^R + 9^(1)^ ]xk = 0, (1) 



— ikt 

with the vacuum-like initial condition Xk — ^75^ in the far past. The coupling to 
the curvature will not be important. Let us seek the solutions of Eq. (|I]) in 
the adiabatic WKB approximation 

a^/\,{t) ^ X,{t) = ^ e-^r-"' + e^^r-^' , (2) 

\/2uj \/2uj 



where the time-dependent frequency is ^^l(t) = ^ + g'^cf)'^ plus a small correction 

~ H'^.H; initially (3k{t) = 0. For \Pk\ ^ 1, an iterative solution 
t . t' 

Pk — I J dt' ^ exp(— 2z / dt"uj{t")) ~ g"^ gives the standard result of the 

—00 —00 

perturbative theory for the particle occupation number nk = |/?fcP- This can 
be interpreted as a pair of infiatons decaying independently into a pair of X' 
particles, {(fxp XX)- The largest total decay rate in the perturbative reheating 
is r < 10~'^^Mp] the upper bound on the perturbative reheating temperature is 
Tr<10^ Gev. 

Recently it was found, however, that in many versions of chaotic infiation 
reheating begins with a stage of parametric resonance^. Indeed, in Eq. (|1]) infia- 
tons act not as separate particles, but as a the coherently oscillating field (p{t). 
The smallness of g^ alone does not necessarily correspond to small The os- 
cillating effective frequency uJk{t) can result in a broad parametric amplification 
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of Xk- In this case the energy is rapidly transferred from the inflaton field to other 
bose fields interacting with it. This process occurs far away from thermal equi- 
librium, and therefore we called it preheating. Because of the non-perturbative 
character of the particle creation, the theory of preheating is rather complicated. 
Reheating is an intermediate stage between inflation and the hot radiation dom- 
inated Universe. Therefore, during or around reheating one can expect initial 
conditions for many parameters of the Big Bang to be settled: the initial temper- 
ature Tr, baryogenesis, relic monopoles and others topological defects, primordial 
black holes, dark matter relics, etc. Recent developments in the theory of reheat- 
ing show that the stage of preheating can significantly change the way we think 
about these problems. 

Among many exciting topics in preheating, I have selected a particular 
issue for these proceedings: how diverse the character of the parametric resonance 
can be in different models of V{(f)), and how it can be unified. This information is 
vital for finding the leading channel of the inflaton decay. The answers are based 
on a study done in collaboration with Andrei Linde (Stanford), Alexei Starobinsky 
(Moscow) and Patrick Greene (Hawaii) This work was supported by NSF 
Grant No. AST95-29-225. 
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2. Stochastic Resonance in ^</>^ -l- \<t>^'}d Theory 

It turns out that the character of the parametric resonance for the general 
equation (|l|) depends on the character of the classical oscillations (pit), i.e., on 
the shape of the effective inflaton potential V(0). We begin with the quadratic 
potential V^(0) = \m^(l?. The background solution corresponds to sinusoidal 
oscillations (j){t) ^ $(t) sin (m^t), with the amplitude $(t) = ■ decreasing 
as the universe expands with the scalar factor a(t) oc t^/^. 

For a the toy model the expansion of the universe is neglected (a(t) = 1, a = 0, 
$ = const), the equation for the fluctuations takes the form of a Mathieu equation 



with a new time variable z = mt and the essential dimensionless coupling param- 

2 ff,2 

eter q = The properties of the solutions of the Mathieu equation are well 

represented by its stability/instability chart; see Fig. 1. There is an exponential 
instability Xk oc expln^z) within the set of resonance bands. This instability 
corresponds to an exponential growth in the occupation number of created par- 
ticles nk{t) oc exp{2fikz). Since ~ 10~^, it is expected that q ~ 10^^ g'^ :§> 1, 
the parametric resonance is broad. The simple standard methods for finding the 
characteristic exponent Hk in the narrow resonance regime g -C 1 do not work 
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Fig. 1. Stability/instability charts of the Mathieu equation in the variables 
(left) and of the Lame equation (|]) in the variables (k^, ^) (right). 
Shaded (unshaded) areas are regions of instability (stability). For the Lame 
resonance, the higher bands shrink into nodes for ^ = Ilill^ll^ when simple 
closed form analytic solutions are found^. The inflaton ^A^"^ self-interaction in 
0{N) theory corresponds to n = 1 for = 1 and n = 3 for N > 1. 



here. However, we made the following observation^: for large q the eigenfunction 
Xk{t) is changing adiabatically between the moments tj, j = 1, 2, 3, where the 
inflaton field is equal to zero (f){tj) = 0. The non-adiabatic changes of Xk{t) occur 
only in the vicinity of tj. Therefore, we expect that the semiclassical solution 
(0) is valid everywhere but around tj. Let the wave Xkit) have the form of the 
WKB solution with the pair of coefficients {a^, j3l) before the scattering at the 
point tj] and the pair {a{'^^ , l3p~^) after the scattering at tj. The interaction term 
around all the points tj is parabolic g'^(f)'^{t) ~ g'^(^'^m?{t — tjY. Therefore, the out- 
going amphtudes (a;^^^, /9^^^) can be expressed through the incoming amplitudes 
(a;^, (51) with help of the well-known reflection Rk and transmission Dk amplitudes 
for scattering from a parabolic potential at tj-. 

( e^"^* cosh A 2e+2< sinh A \ / ai \ , _.,2 

= [ -^e-^-i sinh A e^- cosh X ) [pij ^ A = e". , (4) 

where = argF (^-y^) + ^ (^1 + In the phase accumulated by the moment 
tj is 9l = f dt uj{t), and a scaled momentum k — 






The growth index /Xfc is defined by the formula n^^"^ = ni exp(27ryU-^). IFrom 
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Fig. 2. Left: the characteristic exponent fik of the Mathieu equation @ 
as a function of k"^ (in units of ^/gm^) for q ^ 10^ (dotted curve). Two 
instabihty bands are shown. The sohd curve is derived analytically with a 
simple (improvable) approximation^. Right: the characteristic exponent fij. of 
the mode Eq. (|l]) in an expanding universe for j = 10. fil. changes dramatically 
with j but within the envelope (bold curve) which is obtained from Eq. (|^) by 
taking there sin^ = ±1. Positive and negative occurrences of /x^ for k < tt^^ 
appear in the proportion 3:1. 

(I), one finds the characteristic exponent in the large occupation number limit 

= i- In (^1 + 2e-"''' - 2 sin 9it e'^""' \/l + e--«'^ . (5) 

For the Mathieu Eq. (^ there is a simple solution of the matrix equation (^ 
for an arbitrary j, which allows us to find the resonance band | tan(^fc — ipk)\ < 
e"'^'^^ and to fix the phase 9{o^, cos {6k,tot) = vT+~e^sin {9k — ^k)- The crucial 
expression here is the phase 6k accumulating between two successive zeros of 
0(t), Ok = dt^k^ + g^\t), so that ek-Vk-^^+ s^i^^q + 9-474) ^ 
^^+0(fi;^). The function /i^ derived with the formulas (^ is independent of j and 
is plotted in Fig. 2 (left panel), which is a slice of the Mathieu stability /instability 
chart for given (large) q. That transparently shows the presence of a sequence of 
stability/instability bands as a function of k. Typical half-width of a resonance 
band is fc^ ~ O.lgm^. The effect of parametric resonance crucially depends on 
the interference of the wave functions, i.e., the phase correlation/anticorrelation 
between successive scatterings at the parabolic potentials. The maximal value 
of fi is reached for the positive interference when sin6'tot = — 1 and is equal to 
/i = Mn (^1 + ^ 0.28. The typical value of corresponds to sinOtot = and 
is equal to fi = 2^ In 3 ~ 0.175. For the negative interference when sin 9tot = 1 



6 



there is no resonance. The angle 9 tot is mainly controlled by the phase Q-'j^. 

Let us consider now an expanding universe. At first glance, the resonant 
solution here is more complicated than the Mathieu equation, which for g ^ 1 has 
only a heuristic meaning for our problem. The most important change is the time- 
dependence of the parameter q = g oc 2j is the number of infiaton oscil- 
lations. For the broad resonance case this parameter significantly varies within a 
few infiaton oscillations; hence, the concept of the static stability/instability chart 
and the results for the Mathieu equation cannot be utilized here. Surprisingly, the 
most interesting case, when the parameter q is large and time-varying, can also 
be treated analytically by the method of successive parabolic scatterings. Indeed, 
the matrix mapping @ for the and (31 is also valid in the case of an expanded 
universe. The phase accumulating between two successive zeros of the infiaton 
field is now ^ |^ + O(k^). For large initial values of q, the phase rotation 
between successive scatterings 59k ~ is much larger than tt for all relevant k. 
Therefore, the phase 6tot in formula (|^) can be considered to be random numbers. 
Since there is no memory of the phases, each mapping can be considered as 
independent of the previous ones. As a result the number of particles nl. and the 
exponent /x^ are random variables. The functional form of /i^ for the stochastic 
resonance is different from that for the broad parametric resonance, in one 
has substitute a random phase 9 in the interval [0, 27i] instead of 9tot- An example 
of a stochastic /i^ is plotted in Fig. 2. 

Contrary to the regular resonance, the stochastic resonance is much broader, there 
are no distinguished stability /instability bands, and for certain values of momenta 
the function fil is negative. During the stochastic resonance regime, this function 
changes dramatically with every half period of the infiaton oscillations. Fig. |2] 
shows that it is incorrect to use the structure of the resonance bands of the static 
Mathieu equation for investigation of the stage of stochastic resonance, unless the 
parameter q drops to a value ~ 0(1). 

3. Structure of the Resonance in j(f)^ + y'^'^X^ Theory 

In the theory with the potential V{4>) = |A0'^ it is convenient to use 
the conformal field = acf) and conformal time variable r = \/\(p^ J Then 
(without decay) the amplitude of ip is constant, (p. The background solution is 
given by an elliptic function ^{r) (p cn [t, a(r) = \f^^'T- The Eq. (|l]) 
for quantum fiuctuations Xk can be simplified in this theory. Indeed, the physical 
momentum, p = ^ is redshifted in the same manner as the background field 
amplitude. The conformal transformation of the mode function Xk{t) = a{t)xk{t) 
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allows one to rewrite Eq. (|I]) as 

X'^+(^K' + ^jcn'{r,^)y, = 0, (6) 

where k"^ = The equation for fluctuations does not depend on the expansion 
of the universe and is completely reduced to the similar problem in Minkowski 
space-time. This is a special feature of the conformal theory ^Xcp'^ + ^g^cp'^x^- The 
mode equation (||) belongs to the class of Lame equations. Since the coefficients 
of the Lame equation are time- independent, we can numerically construct 
the stability/instability chart in the variables (^k^,^^, plotted in Fig. 1. The 
combination of parameters g'^/X ultimately defines the structure of the parametric 
resonance in this theory. This means that the condition of a broad parametric 
resonance does not require a large initial amplitude of the infiaton field, as for the 
quadratic potential. The strength of the resonance depends non-monotonically 
on this parameter. 

4. Unified Picture for + ^ ^02^2 Theory 

Preheating in the theory + is efficient only if ^ m. The 

infiaton amplitude $ is extremely large immediately after infiation, ~ lO^^Mp, 
and later it decreases as $ ~ The parameter q, which controls the resonance, 
rapidly changes. As a result, the broad parametric resonance regime in this 
theory is a stochastic process. In the theory + y0^x^ the infiaton field also 
decreases in an expanding universe, but it does not make the resonance stochastic 
because all variables scale in the same way as due to the conformal invariance 

2 

of the model. The parameter which controls the strength of the resonance is y. 
What is the relation between these two theories? Neither of these two theories 
is completely general. In the theory of the massive scalar field one may expect 
terms |0'^ to appear because of radiative corrections. On the other hand, in 
many realistic theories the effective potential is quadratic with respect to near 
the minimum of the effective potential. 

Let us consider the theory ^0^ + f 0^ + ^^^X^ which embraces both limits. 
One may expect that for 3> the mass term ^0^ does not affect the fre- 
quency of oscillations of the infiaton field 0, and the parametric resonance in this 
theory occurs in the same way as in the conformal theory, whereas for <^ the 
resonance develops as in the massive theory. Therefore, a naive thought would 
be that as the amplitude $ decreases from 10~^Mp to the parametric reso- 
nance can be described with the Lame stability/instability chart. It turns out 
that the actual criterion for the regular, stable resonance is that the phase of the 
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Fig. 3. Schematic representation of different regimes wliich are possible in 
the theory + + ^^i various relations between g'^ and A in an 

expanding universe. 



Xfc-wave between two zeros of inflaton oscillations, 6k = J ^ dT\J h? + ^^'^{t) is 
not rotating. The phase rotation within an inflaton half-period is 66k — ^^^^f^ . 
This means that despite the subdominant contribution of the mass term, for large 

2 

^ its presence can modify the nature of the Lame resonance by turning it into 

2 2 

stochastic process if $ < ^ fxlT' Different regimes of parametric resonance are 
shown in Fig. p]. Immediately after inflation the amplitude $ is greater than 
If < ^^'^ , the x-particles are produced in the regular resonance regime 
until $(t) decreases to after which the resonance occurs as in the theory 

+ ^(f^X^ ■ The resonance never becomes stochastic. If ^ > ^^^^ ^ the reso- 
nance originally develops as in the conformal theory, but with a decrease of 
the resonance becomes stochastic. Again, for < ^ the resonance occurs as 

in the theory ^0^ + y^^X^- There is a transient regime of intermittancy between 
regular (stable) and stochastic resonances. In all cases the resonance eventually 
disappears when the field becomes sufficiently small and is replaced by the 
perturbative decay. Perturbative reheating can be complete only if there is sym- 
metry breaking in the theory, or additional interactions with fermions^'*^. 
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